Critical statistics in a power— law random banded matrix ensemble 
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We investigate the statistical properties of the eigenvalues and eigenvectors in a random matrix 
ensemble with Hij ~ |i — It is known that this model shows a localization-delocalization 

transition (LDT) as a function of the parameter fj,. The model is critical at /i = 1 and the eigenstates 
are multifractals. Based on numerical simulations we demonstrate that the spectral statistics at 
criticality differs from semi-Poisson statistics which is expected to be a general feature of systems 
exhibiting a LDT or 'weak chaos'. 



In a recent paper Bogomolny et azB have investigated 
a number of dynamical systems and found remarkable 
similarities between the spectral statistics of pseudointe- 
grable billiards and the critical statistics in the Anderson 
model at the mobility edge. The latter haa^hfien inves- 
tigated numerically by a number of authorsoiScl. At the 
metal-insulator transition (MIT) properties intermedi- 
ate betiwseji those predicted by random matrix theory 
(RMT)EJ^EJ and those of uncorrelated spectra with Pois- 
son statistics were found. On one hand the spectral 
statistics for small differences in energy is reminiscent 
of the universal level repulsion in metals, i.e. the near- 
est neighbor level spacing distribution P(s) behaves as 
P{s) ^ for s ^ 0, with (3 — 1, 2, or 4 for orthogo- 
nal, unitary, or symplectic symmetry of the system. On 
the other hand for s 3> 1, In P(s) ~ —as with a con- 
stant a depending weakly on /? but strongly on the di- 
mensionality d, which is reminiscent of uncorrelated en- 
ergy spectra in Anderson insulators. These results were 
found in numerical simulations on c?-dimensional cubic 
systems using periodic boundary conditions (BC). Re- 
cently Braun et aln discovered that the shape of P{s) 
depends strongly on the choice of the BCs. Upon aver- 
aging over BCs they concluded that for orthogonal sym- 
metry P{s) is very close to the semi-Poissonian form 



P(s) = Ase 



-2s 
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The energy level statistics is universal at the MIT in the 
sense that the level number variance Y?{N) = {{SN)^) 
is proportional to the mean number of levels A'^ ^ 1 and 
the coefficient x is independent of the BCs. Dependence 
of the P(s) function on the BCs in d = 2 quantum-Hall 
systemsQ and appearance of the semi-Poisson statistics 
for a, d = 2 Anderson model with symplectic symmctryu 
obtained as an average over the BCs have also been re- 
ported since then. Interestingly similar spectral statistics 
has been obtained for the case of two interacting particles 
in a one dimensional disordered system at the interaction 
strength producing maximal mixing of the noninteracting 
basis. 

At the transition point (MIT) the statistical proper- 
ties of the spectra and of the eigenstates are linked to 



one another. A remarkable relation between the level 
compressibility x ^^nd the density correlation dimension 
of the eigenstates D2 has been derived int3. 
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The dimension D2 describes how the fourth moment of 
the components of an eisepfunction ip scales with the lin- 
ear length of the systemtj. Generally, the 2p-th moment 
scales as X^idV'il^^} oc L^^p(p~^\ where in the case of 
multifractality Dp is a nonlinear function of p. 

The same parameter D2 describes the scaling of the 
probability overlap of two states with an energy sep- 
aration substantially exceeding the mean level spacing; 
hence the iiame of density correlation dimension. It has 
been showrO that the heavily fluctuating local densities 
still produce a considerable overlap so that at the MIT 
level repulsion is still present. 

In Ref.EJ quantum chaotic systems were numerically 
compared to one of the simplest models that pro- 
vides semi-Poissonian statistics: the short-range plasma 
model (SRPM). The model describes N levels thaLjcpel 
each other logarithmically as in conventional RMTLJ, but 
with the interaction restricted to nearest neighbors only. 
For this SRPM many quantities can be computed ana- 
lytically in the large N limit. The spacing distribution is 
given by Eq. (|^) , the two-level correlation function reads 



P(s) 



The corresponding level number variance, 
I]2(L) 



2 8^ ' 



(3) 



(4) 



leads to the level compressibility x = 1/2 which is also 
the upper limit assumed by the right hand side of Eq. 
(^) . Another way of obtaining semi-Poissonian statistics 
is to leave out every other element from an otherwise un- 
correlated sequence of levelS|r-|Such a 'daisy' model has 
been studied recently in RefJla. 
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By allowing the logarithmic repulsion to act without 
limits in the SRPM one recovers the RMT result in which 
the level spacing distribution is well approximated by the 
Wigner-surmise 

P(,)^|,e-^^\ (5) 

and the the two-level correlation function is given 

dc(s) 



R{s) = 1 - c\s) 



ds 



c{s')ds' 



(6) 



with c(s) = sin(7rs) /ns. From this correlation function 
the level number variance of standard RMT (up to 1/L 
corrections) follows 



ln(27rL) + 1 + 7 



(7) 



where 7 = 0.5772... denotes Euler's constant. The level 
compressibility vanishes; the levels nan be thought of as 
particles of an incompressible fluidEjO. 

In the present paper we investigate the critical spectral 
statistics and the multifractalilty of the eigenstates in a 
random matrix model orisiinally proposed by Mirlin et 
alx3 and later discussed inEZl"E2l. The N x N matrices in 
this model are real symmetric and all entries are drawn 
from a normal distribution with zero mean, (Ha) = 0. 



The variance depends 
ment from the diagona' 



iC distance of the matrix ele- 



m,y)^[l + i\^-3\/By'^]-K (8) 

In Ref.Q it has been shown using field theoretical meth- 
ods that for a fixed _B 3> 1 the statistical properties of 
such matrices for /i < 1 resemble those of RMT. On the 
other hand values /i > 1 in the limit N ^ 00 lead to 
uncorrelated spectra similarly as in the case of banded 
random matriceaSd. The case = 1 was proven to be of 
special importance for it produces critical (multifractal) 
eigenstates and critical statistics. Due to the simplicity 
of the basic model and the possibility of analytical treat- 
ment of a MET further detjails have been revealed recently 
by KravtsovEa and MirlirJlj. 

Here we compare the statistical properties of energy 
spectra and eigenfunction properties of this random ma- 
trix model at fi = 1 with the semi-Poisson statistics of 
the SRPM in a regime so far inaccesable to field theoreti- 
cal methods, at B — 1. We study the shape of the nearest 
neighbor spacing distribution P{s), the two-level correla- 
tion function i?(s) and obtain the spectral comressibility 
X from the asymptotic behavior of the number variance 
E^. We also show the multifractality of the eigenstates 
and provide the correlation dimension i?2. We will show 
that for i? = 1 relation (||) is satisfied provided that the 
spectra and the eigenstates are limited to a small portion 
around the middle of the band. 

In our numerical investigation we have collected the 
spectra of x matrices for N = 800 .. . 4800. The 



power law nature of the problem did not allow the use of 
efficient algorithms (e.g. the Lanczos-algorithm) usually 
applied for the study of the MIT in the Anderson model. 
We have performed a very careful unfolding since the 
N ^ 00 and L — > 00 limits are approached very slowly. 
Based on the values of N we considered, we were able to 
extrapolate the expected behavior of several quanities in 
the N —>■ (XI limit. In all of our subsequent discussion we 
limit ourselves to the middle half of the spectra both for 
the eigenvalues and for the eigenvectors. 




FIG. 1. Density of states of the random matrix ensemble 
at /i = 1 and B = 1 rescaled with a/TV. T he continuous line is 
the semi-circle law n{x) = -y/l — {x/A)^ /2n. The inset shows 
the density of states before scaling. 
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FIG. 2. Level spacing distribution P{s) on a log-log scale. 
The continuous line shows the Wigner-surmize w and the 
dashed line represents the semi-Poisson distribution (|l|) . The 
insert shows the same functions on a lin-log scale. 



First the density of states of the model is plotted in 
Fig. The rescaled function is different from the semi- 
circle law. In Fig. || we show that the P{s) is independent 
of N and shows a non-negligible deviation from the semi- 
Poissonian (|^) form. The deviation is seen both for the 
low-s and the large-s part of the function. The low-s 
part of the P(s) is of the form P(s) ^ s, but with a 
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different prefactor. The large-s behavior differs consid- 
erably from the In P(s) ~ —2s form expected from the 
semi-Poisson distribution (|l|). The deviation from the 
semi-Poisson statistics seems to persist in the N oo 
hmit. The 'peakedness' parametem {q = {s'^)~^) of the 
P(s) converges to a value q — 0.7342 ±0.0009 larger than 
that calculated for the semi-Poisson distribution (|l]) for 
which q = 2/3. We see that the shape of the P{s) is 
indeed intermediate between the semi-Poisson and the 
Wigner-surmise. The spacing distribution in a similas. 
random matrix model has been studied by NishigakiE3 
whose result applied for our case is q — 0.7624 which is 
close to our value. 

Fig. H shows that the level number variance S^(i) con- 
verges to a function that has a linear part, S^(L) ~ xL 
with a slope less than unity, x < Ij showing a nonzero 
compressibility of the spectra. However, this compress- 
ibility seems to converge to a value smaller than the 
X = 0.5 of the SRPM. This behavior is reflected also 
in the two-level correlation function that shows devia- 
tions from Eq. (||) , wWilei it is better described by the 
function calculated irHH. It is of the form of (^) with 
c(s) = 0.25 sin(7rs)/ sinh(7rs/4). Tbs-same kernel applies 
also for models discussed in Refs.E3'Ej. 
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FIG. 3. Level number variance E^(iV) as a function of the 
unfolded N = L. The continuous line shows the RMT case, 
Eq. (0), and the dashed line represents a semi-Poisson behav- 
ior of the form Eq. (^). In the insert we compare the two-level 
correlation function with the one obtained from the SRPM 
Eq. (^) (dashed line). The data are a little-JDit better de- 
scirbed by the function obtained by Kravtsovtj (solid curve). 
The RMT function, Eq. (^), is given as a dashed-dotted curve. 

We have checked the validity of the relation (|^) by com- 
paring the level compressibility obtained from the energy 
spectra with the value extracted from the multifractal di- 
mension D2 of the eigenstates. The results are shown in 
Fig.^. As — > 00 the data extrapolate nicely to a com- 
mon value of X = 0.169 ±0.019 that is less than the value 
0.5 expected from the SRPM and is also smaller than the 
value 0.27 found at the Anderson transitions. Our result 
is in full aggreement with the re peait . analytical estimate 
of X = l/(27r) obtained in Refs.EJIlj and also earlier in 



an analogous random matrix modeES. Note that the D2 
value is obtained as an average of the states in the mid- 
dle of the band. The inset of Fig.^ shows that indeed the 
variance of D2 decreases with increasing A^ while taking 
the full set of eigenstates results in a distribution of these 
exponents that remains broad even in the N 00 limit 
a phenomenon already noticed inEj and discussed incj. 
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FIG. 4. Level compressibility as a function of the loga- 
rithm of inverse matrix size, calculated from the eigenvalue 
spectra (open squares) and from the multifractal proper- 
ties of the eigenstates (crosses) [c.f. Eq. (^)]. The dashed 
horizontal line corresponds to the theoretical estimate of 
X = l/(27r) — 0.159154 .... The errorbars are calculated from 
the variance of the correlation dimension, a{D2), the errorbar 
of the data from the spectrum is smaller than the size of the 
symbol. The insert shows the A^'-dependence of (7(1)2) if the 
averaging is over the full band (circles) or the middle half of 
the band (crosses). 
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FIG. 5. Distribution function of the eigenvector compo- 
nents P{y) where y = {InQ — m)/^/a, m is the mean and a 
is the variance of InQ. The dashed-dotted line stands for the 
power law asymptotics with a power of 1.6. The inset shows 
a close to Gaussian form for P{y)- 



In Table || we collected the relevant quantites obtained 
from extrapolations to the thermodynamical limit. The 
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values of D2 and ao were obtained as an average of the 
corresponding exponents over the ensemble of states in 
the middle half of the spectrum. The D2 and ao for each 
state wece obtained using the standard box counting al- 
gorithmllj. 

A further evidence of the multifractal behavior of the 
eigenstates is also presented in the distribution function 
of the eigenvector components P{Q) with Q = It is 
broad at all length scales, i.e. a close to log-normal form 
is expected. In Fig. ^ the inset shows a very nice nor- 
mal distribution of the variable In Q plotted after rescal- 
ing with the mean and the variance of it. However, the 
figure clearly shows deviation from the log-normal form 
especially for the low-y part of the distribution that may 
be described with a power law tail of the form y^'^. Sim- 
ilar deviations have already been detected and studied 
at the quantum-Hall transitioncj and arCjjdearly seen at 
the Anderson-transition in d = 3 systemsE3, as well. 

In summary we have investigated the spectral prop- 
erties of a random matrix ensemble with entries decay- 
ing away from the diagonal in a power-law fashion. We 
found that the spacing distribution function is similar to 
but deviates significantly from the semi-Poisson distribu- 
tion. We have also compared the numerically obtained 
functions S^(L) and R{s) with the analytically known 
ones of the SRPM which produces semi-Poissonian P{s). 
These functions as well as the extrapolated value of the 
level compressibility x differ from those obtained for the 
SRPM. The latter agrees well with analytical results. We 
also confirm the existense of multifractal states that are 
characterized by a correlation dimension Z?2- According 
to our numerical results the correlation dimension satis- 
fies the relation with x obtained from the spectra. 
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TABLE I. Relevant parameters obtained in the N 00 
limit from the spectral and eigenvector statistics of the states 
in the middle half of the band. 





quantity 


origin 


value 


eigenvalues 


<? 


is'}-' 


0.7342 ± 0.0009 




X 




0.1732 ± 0.0029 


eigenvectors 


D2 


{D2) 


0.6617 ± 0.0389 




X 


(l-52)/2 


0.1691 ±0.0195 
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